Sensing of Fluctuating Nanoscale Magnetic Fields Using NV Centres in Diamond 
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New magnetometry techniques based on Nitrogen- Vacancy (NV) defects in diamond allow for 
the imaging of static (DC) and oscillatory (AC) nanoscopic magnetic systems. However, these 
techniques require accurate knowledge and control of the sample dynamics, and are thus limited 
in their ability to image fields arising from rapidly fluctuating (FC) environments. We show here 
that FC fields place restrictions on the DC field sensitivity of an NV qubit magnetometer, and that 
by probing the dephasing rate of the qubit in a magnetic FC environment, we are able to measure 
fluctuation rates and RMS field strengths that would be otherwise inaccessible with the use of DC 
and AC magnetometry techniques. FC sensitivities are shown to be comparable to those of AC 
fields, whilst requiring no additional experimental overheads or control over the sample. 

PACS numbers: 03.65.Yz, 07.55.Ge, 07.79.-v 



The exploitation of controlled quantum systems as 
ultra-sensitive nanoscale detectors has tremendous po- 
tential to advance our understanding of complex pro- 
cesses occurring in biological and condensed-matter sys- 
tems at molecular and atomic scales 0, 0, 0] . The strin- 
gent requirements for high sensitivity and spatial reso- 
lution has led to suggestions of using spin-based quan- 
tum systems as nanoscale magnetometers Q , or of imag- 
ing through detection of sample induced decoherence [5J . 
One particularly attractive physical platform to imple- 
ment these ideas is the Nitrogen- Vacancy (NV) centre in 
diamond [Fig. [Ha)], chosen for its long coherence times 
at room temperature and convenient optical readout of 
the spin state [fj [Fig.QJb)]. As such, NV centres have 
been the focus of recent proposals to image static (DC) 
and oscillating (AC) magnetic fields @, a|> which have 
since been demonstrated experimentally 2, 3, HI ]. 

However, many important biological and condensed 
matter systems exhibit non-sinusoidal fluctuating mag- 
netic fields with extremely low or zero mean values 
[Fig.QJd)]. An important question is therefore to what 
extent these quantum based magnetometry techniques 
are applicable to such situations. In this paper we ad- 
dress this by quantifying the detection sensitivities for 
these modes for samples with fluctuations characterized 
by the RMS field and dominant spectral frequency. The 
results indicate that by probing the dephasing rate of 
a spin qubit placed in such environments one can char- 
acterize the underlying fluctuation rates and RMS field 
strengths that would be otherwise inaccessible with the 
use of DC and AC magnetometry techniques, thereby 
opening the way for non-invasive nanoscale imaging of a 
range of biological and condensed matter systems. 

The theory behind the detection of magnetic fields us- 
ing quantum systems is heavily reliant on the phase esti- 
mation program of quantum metrology, particularly the 
determination of coupling parameters that are constant 



(a) 



(c) 

2,88 GHz 
microwave 
control 




(b) 




> y 




It)- 




*— IBu) - 


la) J 


f .2.88 GHz 1 



Diamond nanocrystal 
containing NV centre 




(d) Magnetic field at NV 



ii 
in 

IV 



(e) 



Ground state population 



■■ V ■ ,v,-.-- ■ 



FIG. 1: (colour online). Schematic of a scanning NV 
qubit magnetometer/decoherence probe for the detection of 
nanoscale field fluctuations, (a) NV-centre diamond lattice 
defect, (b) NV spin detection through optical excitation and 
emission cycle, (c) Microwave control of the spin state of 
the NV centre and 531 nm optical pulse for read-out. (d) 
Magnetic field signals B(t) at the NV probe corresponding to 
regions I-IV of an inhomogeneous test sample with different 
fluctuation amplitudes and frequency spectra, (e) The corre- 
sponding qubit excited state populations P(t) show that the 
regions can be distinguished by the dephasing information: 
I: Strong, rapid fluctuations — » fast exponential dephasing. 
II: Strong, slow fluctuations — > fast Gaussian dephasing. Ill: 
Weak, rapid fluctuations — > slow exponential dephasing. IV: 
Weak, slow fluctuations — > slow Gaussian dephasing. 



in time. In the context of DC magnetometry, this corre- 
sponds to measurement of the first moment (the mean) 
of the magnetic field strength. For zero mean fields, com- 
plex microwave control pulse sequences are necessary. 
For fields exhibiting oscillatory (AC) time dependence 
with which either a spin-echo or Carr-Purcell-Meiboom- 
Gill [13] sequence may be synchronised, sensitivities are 
predicted to be as low as 3 nTHz -1 / 2 [gj, based on the 
standard quantum limit. Excellent agreement between 
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theory and experiment has been demonstrated in 11 1. 
Such techniques require accurate knowledge of the field 
dynamics which may not be available, or more commonly, 
the field strength may exhibit a stochastic time depen- 
dence. Examples include nuclear dipole fields of ion chan- 
nels [l3| [FigJEfa)] and lipid bi-layers in biological cell 
membranes [1J], Overhauser fields in Ga-As quantum 
dots [HI, and even self-diffusing water molecules 3, 17 1 
[Fig.fJtb)]. i n what follows, we investigate the effects of a 
more general fluctuating (FC) field on the dephasing of a 
spin qubit as the primary detection mechanism, and the 
implications for the characterisation of the magnetic field 
from the surrounding environment. In this sense, we are 
estimating the second moment of the environmental field 
strength, and the corresponding temporal dynamics. 

A spin qubit placed in a randomly fluctuating mag- 
netic environment will experience a complex sequence of 
phase kicks, leading to an eventual dephasing of the pop- 
ulation spectrum. For an NV centre, this will be in ad- 
dition to the intrinsic sources of dephasing, which are 
due to paramagnetic impurities in the diamond lattice 
[l8j |. The dephasing rate can be quantified via repeated 
projective measurements of the qubit state, and the cor- 
responding dephasing envelope, T>(t), can be determined 
via a suitably chosen quantum state reconstruction tech- 
nique. We use the technique of Hamiltonian characterisa- 
tion [l9l | rather than quantum tomography techniques, as 
it requires only a single measurement basis yet is robust 
in the presence of dephasing • 

The motivation for the environment model used here 
comes from consideration of magnetic dipoles in motion. 
Other models in which a two level system is coupled to a 
bath of bistable fluctuators have been previously consid 
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261 ] . These models, however, do 
not capture the dephasing effects due to gradual transi- 
tions between environmental states in slowly fluctuating 
fields. Later we will show this to be of particular impor- 
tance in the case of spin-echo based experiments. Addi- 
tionally, these models require a large number of fluctua- 
tors to model a continuous signal. In contrast, we wish to 
consider the dephasing effects of small numbers of spins 
in motion. 

Consider a qubit with gyromagnetic ratio 7 P undergo- 



ing a f 



Ramsey sequence [27| in the presence of 



a classically fluctuating magnetic field, B{t). An exam- 
ple of a fluctuating magnetic field due to a uni-directional 
spin current [Fig.[2ja)], and that of a bath of self-diffusing 
spins [Fig.[2^b)]. The field has mean (B) = Bq, standard 
deviation y/ (B 2 ) — (B) 2 = B 1 , and typical fluctuation 
rate f e = l/r e , where r e is the characteristic correla- 
tion time of the external field [Fig.[2]Jc)]. This gives rise 
to two natural frequency scales, given by ujq = 7 P -Bn and 
lu' = jpB'. The average precession frequency of the qubit 
is set by wo, and is found to be decoupled from all de- 
phasing effects for cases where u',f e <C ujq. Additional 
relaxation processes may dominate the qubit evolution 
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FIG. 2: (colour online). Typical magnetic field, B(t)/B', for 
(a) a channel of dipoles in unidirectional motion, and (b) a 
self diffusing dipole bath, (c) Temporal correlation function 
(B(t)B(t'))/{B 2 ) . Time axes are rescaled by r e . 



when this condition is violated, however such cases are 
not considered here since we are interested in the char- 
acterisation of weak magnetic fields. The nature of the 
dephasing felt by the qubit will depend on the fluctu- 
ation rate of the environment, / e , or more specifically 
the magnitude of the quantity defined by = f e /^' ■ In 
the case of ^ 1, or fast-fluctuation limit (FFL), the 
qubit will experience many environmental fluctuations 
during its natural timescale. Whilst B' need not nec- 
essarily be normally distributed, the accumulated phase 
error of the qubit at some time t l// e will be nor- 
mally distributed by way of the Central Limit Theorem. 
As such, the variance of the phase error at time t ^> 1 /u>' 
will be (A</> 2 ) ~ t"f 2 B' 2 / / e , giving rise to an FFL dephas- 
ing rate of 
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This is akin to the motional narrowing result from NMR 
[28| and reproduces the ubiquitous exponential dephasing 
envelope given by V fast (t) = exp {-T iast t). 

In the slow-fluctuation limit (SFL), where 8 < 1, 
we note that the magnetic field may be locally approxi- 
mated by a Taylor expansion in t about some initial time 
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where each of the ak has a specific statistical distribution 
containing information about the fcth order derivative of 
B(t), and thus gives rise to a different dephasing channel. 

For the special case where the au are normally dis- 
tributed with mean and variance a 2 (as consistent 
with random dipole motion), the resulting density ma- 
trix following the free evolution time r, but prior to the 
second 7r/2 pulse is defined by p n = p 2 2 = 1/2, and 



dow( r ); where 
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Thus we see the emergence of a hierarchy of dephasing 
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and beating channels, with the dephasing rates and beat 
frequencies of the fc th channel given by 
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respectively. In the case of the zeroth order channel this 
corresponds to the rigid lattice result from NMR [28| , and 
we have Cq = (B 2 ) — (B) 2 . This effect will be suppressed 
by a spin echo pulse sequence. For the first order channel, 
we may approximate u\ ~ ((B 2 ) — (B) 2 ) f 2 . 

The relative contributions of each channel to the over- 
all dephasing rate of the qubit depend explicitly on 
the dynamics of the field, however, it should be noted 
that dominance of the zeroth order channel (ie T?^ > 

^siow» y 3 — 1) ^ s a necessary and sufficient condition 
for the system to exist in the slow fluctuation regime, 
9 « 1. This justifies the use of the Taylor expansion, 
since the resulting polynomial may be well approximated 
by a low-order truncation. 

The intermediate regime of ~ 1 is more compli- 
cated. Fig.[3la) shows dephasing envelopes for various 
values of 0. For times much longer that r e , pure expo- 
nential dephasing behaviour is observed in all cases (with 
dephasing rate Ff ast ), however fast fluctuating environ- 
ments still exhibit slow (Gaussian) dephasing behaviour 
on timescales r where lu't < >/2/Q. If is large, con- 

(k) 

tributions to T> from the Fg low will decay rapidly. The 
abrupt transition from X> 8 i ow — * £>f as t is shown more 
clearly in the corresponding insert. 

For the purpose of comparison with existing spin- 
based magnetometer proposals, we take the NV centre 
as our example qubit. The Hamiltonian used to de- 
scribe the time evolution of an NV-centre is given by 
H = S • D • S + fi/7 p B • S + 7i thcr, where 7i ther describes 
higher order effects such as hyper-fine splitting, interac- 
tion with optical fields, etc. which can be ignored in the 
present context. We consider weak external fields such 
that O (tvypB ■ S) < O (S • D • S), thereby ensuring the 
crystal-field splitting tensor, D, sets the quantisation axis 
of the NV centre, and that u>' *C loq (even in the FFL). 

The shot-noise-limited DC magnetic field sensitivity 
for an NV-based magnetometer subject to a Ramsey- 
style pulse sequence is given by [§] 



rjdc 



T 



1 



(6) 



where yT and C represent the combined effects of spin 
projection and photon shot noise for N s measurements 
[C — ► 1 for the ideal case), r is the free evolution time 
of the qubit in a given experiment, and T = N s t is the 
total averaging time for N s such experiments. Dephas- 
ing times due to the interaction of the NV centre with 




FIG. 3: (colour online), (a) Plot of simulated dephasing en- 
velopes for N s — 10 4 runs, showing agreement with Eqs. Q] 
& [2] Time is in units of {p/ v B') 1 . (insert) Zoomed plot 
showing that fast fluctuating environments still exhibit non- 
exponential dephasing for short timescales r : ui't < \/2/Q. 
(b) log- log plot of DC magnetic field sensitivity, r/dc, as a 
function of / e for different contours of B' . Assumed parame- 
ter values are T 2 * = 1 fis, andC = 0.3. 



nearby paramagnetic lattice impurities will in general be 
different for different centres and will thus require indi- 
vidual characterisation. For comparison with [8(], we use 
the commonly accepted value of t = T 2 * ~ 1 fis. 

We emphasise here that Eqn. [6] applies solely to the 
imaging of DC magnetic fields where the dephasing of 
the qubit is exclusively due to intrinsic crystal effects. 
If the sample being imaged produces a fluctuating field 
of sufficient amplitude, the dephasing time (1/r) may 
be shorter than T 2 *, resulting in poorer static field sensi- 
tivity. In this context, ?7dc refers to the sensitivity with 
which the mean field, (B), may be measured as the field 
fluctuates over the course of the experiment. To gain 
insight into the effect of fluctuating magnetic fields on 
the DC field sensitivity, we consider again a | — r — | 
sequence. The DC sensitivity as a function of B' and f e 
is shown in Fig. 0(b). From this, we see that fluctuating 
environments can have a dramatic effect on the DC field 
sensitivity of an NV based magnetometer, depending on 
both field strength and fluctuation frequency. 

We now turn our attention to the the sensing of mag- 
netic field fluctuations themselves, including the case of 
zero-mean magnetic fields. Using coherent control tech- 
niques (spin-echo for example), we may extend the de- 
phasing time of the NV-centre to ~ 300 /is, as dictated 
by the 1.1% carbon-13 content in the lattice. The case 
of perfectly oscillatory magnetic fields has been consid- 
ered in detail in [8j, in which AC sensitivities may be as 
low as 3 nT Hz~ 1/2 [Fig-Efb)]. This technique requires 
the 7T pulse to coincide with the first zero-crossing of the 
magnetic field, which requires an accurate knowledge of 
the oscillation frequency. This may prove difficult unless 
the frequency is externally controlled. Furthermore, ac- 
curate control will become difficult at high frequencies, 
and this may lead to further dephasing. 

Rather than considering an AC field, we now study 
the magnetometer's sensitivity to a more general fluc- 
tuating field via consideration of the induced dephasing 
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rate [5(. For a f - £ — it — \ — | pulse sequence, the 
probe will show decreased sensitivity to environments ex- 
hibiting fluctuation frequencies, f e less than 1/r. For 
fast fluctuating fields, the effect will be negligible. For 
the imaging of slowly fluctuating fields, this may appear 

problematic, however complete insensitivity only comes 

(k) 

with f e — ► 0. A spin echo sequence will modify the £>g low 

(fc) , k A _ 2 -fc 



4 (a) Fluctuating field source 
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^sSw tnus on ^y * ne effects of the 
zcroth order dephasing channel will vanish. 

Perturbations on the dephasing rate may be measured 
from (1 — T>) min — CXP q T J^ ^ i where n describes the 
shape of the spin-echo dephasing envelope as dictated 
by the presence of carbon-13 nuclei in the lattice, which 
for present purposes may be taken as n — 3 [18]. This 
implies an optimal free-evolution time of r ~ Tij \f§. 

Thus we find that perturbations on the 1/T 2 dephas- 
ing rate as slow as 200 Hz for exponential dephasing and 
800 Hz for Gaussian dephasing may be detected by this 
method after 1 s of averaging time. By performing mea- 
surements of the total dephasing rate, T, both the field 
variance and average fluctuation rate may be inferred 
from equations [T]&[4j Of course, the question remains of 
which fluctuation regime a given sample system resides 
in. In the absence of any prior knowledge of the environ- 
ment being measured, this question may be answered via 
determination of the shape of the dephasing envelope, a 
task to which the Hamiltonian Characterisation method 
is well suited 

The optimal fluctuating magnetic field sensitivity will 
occur when the ~ 1 condition is satisfied, since this en- 
sures maximal dephasing for a given field variance. Con- 
sidering the special case of pure exponential decay, we 
therefore expect an optimal sensitivity of r]f c — c ^ , 

i/2 P 

which for C = 0.3 gives r/f c = 1.7nTHz ' . In practice 
however, such sensitivity may be difficult to realise due 
to memory effects in the fluctuating environment. For 
systems that satisfy ^> 1, thus exhibiting long-time 
exponential dephasing behaviour, Gaussian dephasing is 
still exhibited for t < l// e [Fig.[3][a)]. For spin-echo ex- 
periments, the effect is worsened as the dominant contri- 
bution to ©slow comes from k = 1. Taking this into con- 
sideration, the minimum resolvable field obtained after 
T = 1 s averaging time is plotted in Fig.[4ja) as a function 
of environment fluctuation frequency. We see that fluctu- 
ating field strengths as low as 4.5 nT may be achievable 
after T = 1 s averaging time, or some N s ~ 3000 mea- 
surements. It is also evident that the qubit will be sen- 
sitive to FC fields fluctuating on timescales much slower 
than I/T2. This is in direct contrast with the AC case, 
which shows poor sensitivity to fields oscillating with pe- 
riods less than T 2 [Fig.g^b)]. 

We have theoretically investigated the effects of a fluc- 
tuating magnetic field on an NV centre spin qubit. This 
analysis was used to place new limits on the sensitivity 
with which the mean field strength may be measured. 
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FIG. 4: (colour online), (a) Minimum resolvable FC field 
strength, B m i u , vs environmental fluctuation rate, / e , for 
T — 1 s averaging time. In contrast to the AC case, an FC de- 
tection requires no prior knowledge of fluctuation timescales. 
(b) Minimum resolvable AC field amplitude, B^ in , vs field 
oscillation frequency, i/ e , for T = 1 s averaging time. Here 
we have assumed that the AC field is initialised in phase with 
the probe qubit, and that the oscillation frequency is known 
accurately enough that the ir pulse of a spin-echo sequence 
coincides with the first zero-crossing of the field. This will 
become increasingly difficult with increasing v e , leading to 
additional sources of dephasing. 



Furthermore, we have built upon the idea of decoherence 
microscopy [5] to theoretically demonstrate the ability 
of an NV centre to measure field strengths and fluctua- 
tion rates of randomly fluctuating magnetic fields. This 
analysis shows that the methods presented here require 
no experimental resources beyond those of existing tech- 
niques, no prior control or knowledge of the external field, 
and thus may be implemented with current technology. 
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